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Abstract 

Non-BPS non-Abelian vortices with CP 1 internal moduli space are 
studied in an N = 2 supersymmetric £7(1) x SU(2) gauge theory with 
adjoint mass terms. For generic internal orientations the classical force 
between two vortices can be attractive or repulsive. On the other 
hand, the mass of the scalars in the theory is always less than that of 
the vector bosons; also, the force between two vortices with the same 
CP 1 orientation is always attractive: for these reasons we interpret 
our model as a non-Abelian generalization of type I superconductors. 
We compute the effective potential in the limit of two well separated 
vortices. It is a function of the distance and of the relative colour- 
flavour orientation of the two vortices; in this limit we find an effective 
description in terms of two interacting CP 1 sigma models. In the limit 
of two coincident vortices we find two different solutions with the same 
topological winding and, for generic values of the parameters, different 
tensions. One of the two solutions is described by a CP 1 effective 
sigma model, while the other is just an Abelian vortex without internal 
degrees of freedom. For generic values of the parameters, one of the 
two solutions is metastable, while there are evidences that the other 
one is truly stable. 



1 Introduction 



According to the ideas of Mandelstam and 't Hooft pQ, confinement in non- 
Abelian gauge theories is due to a dual Meissner effect. The electric flux 
between two heavy electric sources is confined by a dual Abrikosov-Nielsen- 
Olesen vortex [2]; the vortex has a constant energy per unit length (tension). 
This leads to a linear potential between the probe charge and anti-charge. 
Due to the difficulties in analyzing strongly interacting non-Abelian gauge 
theories, this picture remained just a nice qualitative scheme for years, which 
could not be justified from first principles. 

A breakthrough in this context was the Seiberg-Witten solution [5] of 
M = 2 super Yang-Mills theories; they found massless monopoles at strong 
coupling. Adding a small M = 2 breaking mass term for the adjoint field, 
the monopoles condense creating dual vortex strings which carry a chromo- 
electric flux. The details of confinement in the Seiberg-Witten scenario are 
indeed quite different from QCD. The SU(N C ) gauge symmetry is sponta- 
neously broken broken to U(l) Nc ~ 1 by the expectation value of an adjoint 
field and the strings of the theory carry an Abelian [/(l)^ -1 charge. A 
careful examination shows that the "hadronic" spectrum is much richer than 
that of QCD [1]. 

Thus it is interesting to study some alternatives to the Abelian Meissner 
effect, with the aim to find some close relatives of QCD. The non-Abelian 
vortex discussed in Refs. [U [6] is an interesting possibility in this direction. 
This solitonic object has first been studied in an Af = 2 U (N c ) gauge theory 
with Nf = N c = N quark hypermultiplets and with a Fayet-Iliopolous term 
in order to keep the theory in the Higgs phase. The squark fields condense 
and break the gauge symmetry; on the other hand the colour-flavour locked 
global symmetry is unbroken in the vacuum. The theory has vortex solitons 
which spontaneously break this SU(N) symmetry to SU (N — 1) x U(l); due 
to the zero modes corresponding to these broken symmetries, the moduli 
space is given by the quotient: 

cp»- su ^ 



SU(N-l) x [/(!)' 



The classical moduli coordinate can be promoted to a field living on the 
vortex worldvolume; in this way vortex solitons in a 3 + f dimensional theory 
can be directly connected with a CP^ -1 sigma model in f + f dimension, 
which describes the macroscopic physics of the flux tube. Some reviews on 



f 



this subject can be found in Ref. [7]; recent developments involve the Seiberg 
duality [8J, the Goddard-Nuits-Olive duality [H], generalizations to SO(N) 
theories [10J, and possible phenomenological applications to high temperature 
Yang-Mills and dense matter [11] . 

The non-Abelian vortex can be studied in many different theoretical set- 
tings; indeed it is possible to start with an M = 1 theory [T2J [T3] or even 
with a non-supersymmetric theory [13]. The details of the effective f + f 
dimensional sigma model are different due to different number of fermions 
and various amount of supersymmetry. Also the number of quantum vacua 
is different, for example there are N vacua in the M = 2 case [151 EH E] an d 
just one vacuum in the non-supersymmetric case [Ti] . 

It is also interesting to study the non-Abelian vortex for higher winding 
numbers. In the M = 2 case, the vortex is a BPS object with a big moduli 
space; as discussed in Ref. [3 [18], in the topological sector with winding k, 
the dimension of the moduli space is 2kN. Some of these moduli correspond 
to the relative and global positions of the component k = 1 vortices; others 
to the global and relative orientations in the internal space. The vortex 
solution and the moduli space for higher winding numbers has been discussed 
in Refs. pH M, M, I2D 122] • 

In this paper, we study the impact on the vortices of the Af = 2 model 
of some mass terms f]o,f]3 for the adjoint fields, which break the extended 
supersymmetry. For concreteness, we will discuss the case N c = Nf = 2. The 
vortex with winding number one is not anymore a BPS object, but still has 
a CP moduli space. On the other hand the physics for vortices with higher 
winding numbers is very different: almost all the flat directions in the moduli 
space are lifted by the parameters r)o, 773. The force between two vortices is 
not as simple as in an Abelian superconductor, where we have attraction for 
type I vortices and repulsion for type II ones [23, 24J. There is a non-trivial 
dependence on the orientations of the two vortices in the internal space. 

Even if the force between two vortices in our model is not attractive for 
all values of the vortex orientations Hi, H2, we have a close resemblance with 
type I Abelian vortices: we find that the scalars of the theory are lighter 
than the vector bosons. Hence if we consider two well separated vortices, we 
have that the prevailing part of the interaction is mediated by scalars and 
not by vectors. Moreover, for n\ = 712 the force is always attractive. We have 
also found evidences that the configurations which minimize the energy are 
always given by two coincident vortices, just as in the type I Abelian case. 
For these reasons we call these objects "non-Abelian type I vortices". 
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In Sect. 2 we describe the theoretical set-up. In Sect. 3 we discuss the 
non-BPS solutions for the vortices which live in the Abelian subset of the 
theory, with emphasis on the sectors which have topological winding 1 and 2. 
In Sect. 4 we study a more general configuration of two coincident vortices, 
both in the BPS case and in the non-BPS potential for the vortex 

moduli space is found for 770,3 7^ 0. In Sect. 5 the interactions between two 
vortices with a large separation distance are studied and the effective vortex 
potential is computed in this limit. In Sect. 6 the worldsheet description of 
the macroscopic physics is discussed, both for a single vortex and for two 
vortices at large distance. In Sect. 7 we conclude the paper and make a 
general discussion. Some aspects specific to the large r/ 3 limit are discussed 
in Appendix A. In Appendix B the BPS equations for two coincident vortices 
for rjo = rj3 = are provided. 

After this work was finished, two papers about interactions of global non- 
Abelian vortices appeared [25J; however the details of these two models are 
quite different from the setting studied in this paper. 

2 Theoretical Set-Up 
2.1 Lagrangian 

For the 17(1) gauge field A° and the SU{2) gauge field A k (k = 1, 2, 3) the 
following conventions are used: 

A, = + \a% (1) 

V, = d, - zy^ - l -A% D,a k = «V + e klm 4a m . 
The field strength is: 

F, v = d,A u - d v A, - % - [Ajy, A k v r k ] , (2) 

which in components (with the convention = F k u r k /2) reads: 

^/L = d^AI — d u A 1 ^ + e^ k A^A k . 



3 



We consider an Af = 2 supersymmetric U(l) x 577(2) gauge theory with 
Afy = 2 hypermultiplets together with the following superpotential: 



V2 



Q f {a + a k r k + v / 2m / )Q / + W (a) + W 3 (a k r k )] , (3) 



where the terms Wq j3 are of the form 

W = -£a + f] Q a 2 , W 3 {a k r k ) = r] 3 a k a k . (4) 

Here we have introduced two real positive mass parameters r] and rj 3 for 
the adjoint scalars which break M = 2 SUSY to M = 1. my is the mass of 
the hypermultiplets Qf,Qf (f = 1,2) and £ is the FI F-term parametei0. 
Vortices in Abelian versions of this theoretical setting have been discussed in 
Refs. [28], [29]. 

This kind of potential naturally arises from the M = 2, SU(3) SQCD 
softly broken with a mass term of the form W = r]TiA 2 . Indeed, when 
the bare masses of the squarks are tuned to special values, there exist true 
quantum vacua in which the non-Abelian gauge symmetry SU(2) x U(l) is 
preserved [6]. The low energy effective theory in these vacua is exactly the 
theory we are studying herej 

The bosonic part of the Lagrangian in Euclidean notation is (we use the 
same symbols for the scalars clS 0X6 used for the corresponding superfields): 



C = / (fx 



+ Tr(V M Q) t (V M g)+Tr(V M Q)(V M g t )+y(Q,g,a fe ,a)], (5) 

where e is the J7(l) gauge coupling and e 3 is SU(2) gauge coupling. The 
potential V is the sum of the following D and F terms: 



V = f (^j jk &a k + ^(Q*T i Q)-Tr&&)\ 



1 A very similar Lagrangian was discussed in Refs. [12], [13]; in that case the FI was 
in the D-term and not in the superpotential. This leads to different physics, the vortex 
is still classically BPS saturated. For a discussion of the different settings that give BPS 
vortices, see Refs. [26]. 

2 These quantum vacua exist only if we have a sufficient number of flavours. In this 
case semilocal vortices may be relevant |27j . 
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+ |(Tr(QtQ)_ T r(QQt 



+ 



+ 



Tr (QQ) - £ + 2r] a 



i 2 

+ -^|(a + rV + v / 2m / )g / | 2 + |(a + rV + v^m^Q}! 2 . (6) 
/=i 

The squark multiplets are kept massless in the remainder of the paper, 

rrif = 0. 

The vacuum of the theory which we are interested in is not changed by the 
parameters 770,3: 



Q = Q = \ll( I ?)» a = > ab = - 

For rjo 7^ 0, the theory has also another classical vacuum: 

Q = Q = 0, a = a k = 0, 
2r/ 



(7) 



(8) 



which "runs away" at infinity for r/ = 0. In what follows, we consider the 
vacuum (j7|) and the non-Abelian vortices therein. 



2.2 Spectrum of the Theory 

The masses of the gauge bosons can easily be read of the Lagrangian: 

M£ (1) =£eg, M 2 SU{2) = ie\. 

The masses of the scalars are given by the eigenvalues Mf of the mass matrix 
(calculated in the vacuum (j7|)): 

where we denote by s& (k = 1, 2, • ■ ■ , 24) the real scalar fields of the theory. 
The calculation is a bit tedious but quite straightforward (a very similar sit- 
uation is discussed in Ref. [12], in the case of a D-term Fayet-Iliopoulos). 
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First of all, there are four zero eigenvalues, which correspond to the scalar 
particles eaten by the Higgs mechanism. There is one real scalar with mass 
M50 = Mu(i) which is in the same M = 1 multiplet as the U(l) massive 
photon and moreover there are also three scalars with a mass Mto = Msu{2) 
in the same multiplet as the non-Abelian vector field. The other mass eigen- 
values are given by 

M 2 S1>S2 = £ejj + e 4 oV 2 ± ^gejj + efog ; (10) 

where the upper sign is for M 51jT1 and the lower sign is for M S2 ,T2- Msi and 
M52 have multiplicity 2; M?\ and Mt% have multiplicity 6. The mass of the 
fermions is obviously the same as the mass of the bosonic degrees of freedom, 
because of the unbroken M = 1 supersymmetry. 

Note that for rjo = 773 = (which is also discussed in Ref. [30]), the mass 
degeneracy of the spectrum is bigger (the particles fit into M = 2 hypermul- 
tiplets). The parameter 770 affects only the masses of the particles which are 
in the same Af = 2 hypermultiplets as the U(l) vector field; 773 affects the 
mass of the particles which are in the non-Abelian vector hypermultiplet. 
The case of M = 2 SQED was studied in Ref. [28]; the results are very similar 
to our £7(1) subsector. 

In the limit r]oe > v^, we find M| x w 2e^g and M| 2 w £ 2 /( 4? 7o)- 
In a similar way, if 77363 3> \/£ the masses become M|i w 2efr?f and 
M|, 2 w (4773 ) . The particles with masses Mg\ and Mti (which in this 
limit correspond to the fields a and a k ) become very massive and decouple 
from the low energy physics. 

For r] eo a/£ an d 7/363 <€i nonetheless we find 

M ks2 = &l ± y/£vo4, M 2 T1)T2 = ie\ ± s/lrfr4. (11) 

Some of the scalars become slightly heavier and some slightly lighter. 

The mass eigenvectors take a quite complicated form for small 770,3, with 
at non-trivial mixing between Q, Q and a, a k . On the contrary they are quite 
simple for large 170,3, because the fields a, a k decouple from the low energy 
physics. The effective Lagrangian for large 770,3 is discussed in Appendix A. 
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3 (p, k) Coincident Vortices 



3.1 Second order equations 

In this section, we will study some special solutions representing coincident 
vortices that live in an Abelian subgroup of the fields of the theory. These 
solutions are parameterized by two positive integers (p, k) ; the topological Z 
winding number is given by w = p + k. This kind of solution gives us the 
most general vortex with winding w = 1 up to a colour-flavour rotation^. 

Due to the symmetry between Q and Q\ for the vortex solution we 
can consistently set Q = @. With this assumption the Euler-Lagrange 
equations of the theory are: 

d^ v + e klm A ¥ F^ = e 2 3 Tr (iQ^r k (V v Q) - i{V v Q)^r k Q) - 

- e k i m {{D v a)\a m + {D v a)ia\ n 



mil 



V%Q = -^r, d^a = -^, D^D,a k = -^r. (12) 

We make the following axial symmetric ansatz (which in the BPS limit re- 
duces to the one of Ref . [6] ) : 



Q 



A i = "^[(P - *) - Ur% A* = - 6 -^i[(p + k)- / (r)], 



ao = A (r), a 3 = A 3 (r), a x = a 2 = 0. (13) 

Notice that the adjoint fields a, a k are non-trivial in the non-BPS model, 
whereas they vanish everywhere when 770,3 are zero (BPS). 

The vacuum of the theory is invariant under the following global colour- 
flavour locked rotations (U G SU{2)c+f)'- 

Q^UQU\ Q^U^QU, a k T k ^U(a k T k )U\ F k u r k - U(F k u r k )Ul (14) 



3 As we will discuss in the next section, for higher winding, we know that it is not the 
most general solutions, at least in the BPS case [5l fT9 l [TB I [20 l [2Tj. 

4 This can be checked using the variables: Q S = (Q + <9 f )/2, Qd = {Q ~ Q f )/2. With 

= 0. 



this variables we easily see that -§q^ 



=0 
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Let us introduce the S 2 coordinate n k , with k = 1, 2, 3 and \n\ = 1: 

n k T k = Ur 3 Ul (15) 

Using the parametrization introduced in Eq. (115)) we can write down the 
expression for a w = 1 vortex with generic orientation n k : 

A° = -^[l-/o], ^ = -^[1-/3]"*, a» = A , a* = n %, 

Q = Q f = *4±^1 + (16) 

It is easy to see that the (1,0) vortex partially break the symmetry in 
Eq. (TH1) ; as a consequence, this object has some internal zero modes associ- 
ated to this breaking. In fact, the vortex leaves a U (1) subgroup of SU {2)c+f 
unbroken, so that zero modes parameterize a CP 1 = SU(2)/U(1) = S 2 . 
The energy with respect to 01,2, /o,3 an d Ao,3 is expressed as 

// f/2 f/2 \/2 \/2 

r *'(^ + ^ + ^ + l + 2W + * ?)+ 

w + mi + /i) + 2/3/0(0? - 01) + * A + A _ { + J+ 



2r 2 2 



+f (0i - ^2 + 2%A 3 ) 2 + ((A + \,)<PiY + ((A - A 3 )0 2 ) 2 J • (17) 

We have to minimize this expression with the appropriate boundary condi- 
tions for each (p, k): 

f 3 (P)=p-k, f o (0)=p + k, / 3 (oo) = 0, / (oo) = 0. 

<Moo) = l, 02(oo) = 1, A (oo) = 0, A 3 (oo) = 0. (18) 

We also find for small r: 

lt x0(rP), 2 ocO(r fc ), A a 0(1), A 3 a 0(1). (19) 

The Euler-Lagrange equations obtained are: 

^-4 = -(/3(0i-0i) + /o(0 2 + ^)), 
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^-4= +«£)+/o(#-^))> 




0i ((A + A 3 ) 2 + e 2 (0 2 + 01 - ^ + S^Ao) + e 2 (0 2 - 0| + 2 m X 3 )) 




02 ((A - A3) 2 + e 2 (0 2 + 2 - £ + 2?? A ) - e|(0 2 - 2 + 2/7 3 A 3 )) 





Aq = e 2 ((g + a 3 )0 2 + (o - a 3 )0 2 + ejfoo(0 2 + 2 - £ + 2r7 A )) 
r 2 

A^ = e| ((fl + a 3 )0 2 - (op - a 3 )02 + ej^C^i ~ 4>l + 2t7 3 A 3 )) 
r 2 



(20) 



It is easy to check that these equations can be obtained substituting the 
ansatz (|T3|) in Eqs. (fl2"|) . This shows that the ansatz is consistent. 

In the following, we will concentrate our effort on the study of the sectors 
with topological winding 1 and 2; in other words, we will discuss the (1,0), 
the (1, 1) and the (2, 0) vortices. In the BPS limit (r] = r] 3 = 0) the tension 
is proportional to the topological winding number (T(i )0 ) = 27r£, = 
T(2,o) — 47r£). For non-BPS solutions, 770,3 7^ 0, we find that the tension is 
always less than the BPS limit. This is because the non-BPS terms in the 
tension formula Eqs. (ITT)) do not give any contribution if we put the BPS 
solutions into the expression (Ao,3 are identically zero for the BPS solutions). 
The non-BPS solutions will of course be a true minimum or saddle point of 
the energy functional, so that their energy will be smaller than that of the 
BPS configurations^. 

For fixed £, the tension of the (1, 1) vortex is a function of only eo,?7o, 
because for this vortex ^3 = and 0i = 02- This is clearly explained by 
the fact that the (1,1) vortex is completely Abelian. On the other hand, the 
tension of the (1, 0) and of the (2, 0) vortex is a non-trivial function of all the 
parameters eo,3, 770,3- 

If we take e 3 = eo and 773 = 770 the vortex becomes easier to study. In this 
case we can use a more convenient basis for the gauge field, which is just the 

5 This also has a clear resemblance to the Abelian case, where type I vortices have a 
smaller energy with respect to the BPS case. 
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sum and the difference of A ^ and A z . The potential V also factorize, and 
takes the form V = Vi(4>%) + V 2 (0 2 )- Each diagonal component of Q does not 
interact with the other ones, and can be treated as an Abelian vortex. For 
the (1,0) vortex we can use the simple ansat20 



while for the (1, 1) and for the (2, 0) vortices we can use 

n-(^y m ^ O-i^ ^ (22) 
Q ~ V 0(r 2 )e^ j ' W " V ) ' 1 J 

The system reduces to the Abelian vortex studied in Ref. [25] • The tension 
of the (1, 1) vortex is exactly twice the tension of the (1, 0) one. In each of the 
U{1) factors, we have type I superconductivity. Since the two U (1) subgroups 
are decoupled, the (1, 0) and (0, 1) vortices do not interact. Furthermore, the 
tension of the (2, 0) vortex is less than twice the tension of the (1, 0) vortex. 

3.2 Numerical Solutions 

At generic eo,3,??o,3 Eqs. ( 120]) have been solved numerically. It is a little 
subtle to solve this system of ordinary differential equations directly. The 
difficulties basically arise because there are many equations; there are also 
subtleties in defining the boundary conditions at oo, because, in general, the 
fields which appear in our ansatz do not correspond to mass eigenstates. In 
order to perform the numerics we found that the method of relaxation is 
very effective. We add an auxiliary time dependence to the profile functions 
w = (foi hi 0i) 4 ) 2, Ao, A3) . At t = we start with some arbitrary functions 
Uj(r, Ojj, the evolution in t is then given by: 

sf = < 23 > 

If the solution converges with time to a static configuration, then at final time 
we have obtained a solution of the equations Ej = 0, which are equivalent to 
Eqs. ([20]). The results for (p, k) = (1, 0), (2, 0), (1, 1) are shown in Fig. GQ 



6 Notice that we cannot impose <f>2 — Vs/2 even for (1, 0) vortex in the generic models. 
7 The choice of the initial conditions is crucial to find convergence. 
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Figure 1: Profile functions /o (solid black), fs (solid blue), 4>i (long dashes, black), 02 
(long dashes, blue), Ao (short dashes, black), A3 (short dashed, blue) for the numerical 
values £ = 2, eo = 1/4, e3 = 1/2, 770 = 773 = 1. In the left panel are shown the profiles for 
the (1, 0) vortex, in the middle, the ones for the (2, 0) and in the right, the ones for the 
(1,1). Note that in this last case /a = A3 = and <\>\ = (j>2- 

It is interesting to compare numerical result for the tension 2T 10 , T 2j0 and 
T 1A of the 2 x (1,0), (2,0) and (1, 1) vortices, respectively (see Fig. Ej). We 
have always found that T 2 o < 2T 10 . This is consistent with the fact that 
at large separation distance, the force between two vortices with the same 
colour-flavour orientation is always attractive (we will discuss this aspect in 
Sect. 4). As can be checked in Fig. [2], three different regimes have been found 
for 71,1: T M < T 2i0 < 2T 1>0 or T 2i0 < T M < 2T 1>0 or T 2j0 < 2T 1>0 < T lyV 




Figure 2: Ti,i (long dashes), T2.0 (short dashes) and 2Ti, (solid) for different values 
of < lu < vt/2, where 773 = r)smui,r)o = rjcosuj. The tension of the BPS 2-vortex is 
normalized to Tbps = 2. In the left figure the numerical values eo = = 1/2, rj = 4, 
£ = 2 are used; in the right figure eo = 1/2, e3 = 1/4, 77 = 4, £ = 2). 

If 770 = 0, 773 7^ 0, the tension of the (1, 1) vortex is found to be the same 
as in the BPS case. The tensions of the 2 x (1,0) and of the (2,0) vortex 
are strictly less than that of the BPS vortices (773 = 0). Hence in this case 
T 2t o < 2T 10 < T11. On the contrary, if 773 = 0,770 7^ 0, we have found 
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Ti i < ^2,0 < 2T 10 for all the numerical values of the couplings that we have 
investigated. 



4 Generic Coincident Vortices 
4.1 The BPS case 

The number of dimensions of the fc-vortex moduli space in U(N C ) M = 2 
gauge theory with Nf = N c = N hypermultiplets has been computed in 
Ref. [5]. The calculation uses the index theorem and the result is 2kN. Thus 
for rjo = i]3 = 0, the moduli space of a 2-vortex configuration is a manifold 
with eight real dimensions. Two of these dimensions are associated with 
the global position of the system; other 2 coordinates are associated with 
the relative position R of the two elementary vortices. The remaining 4 
coordinates are associated with the orientation of the system in the colour- 
flavour space. In this section, we will write an ansatz for the case of coincident 
vortices (R — 0), and we will show that it is non-trivially consistent with the 
second order equations of the theory. We will correct a technical mistake in 
Ref. [20J, where the problem was studied using first order BPS equations. 
Modulo a global SU{2) rotation we can parameterize a subset of the 



4) = -^(2 - fo), 4) = -^((i + cos «) - /»). 

4 ) = -^(sina)(cos^)(l-^), Afa = +^(sine*)(sin ¥ >)(l - g). (25) 



We have introduced here the profile functions /«i(r), K%(r), /^(r), K^(r) for the 
squark scalars and /o( r )> /s( r ), d( r ) f° r the gauge field. For r — > oo all the 

8 The ansatz for Q used in Ref. [20] is: K\ = kz± z%, K2 = kzi, = KZ2, ka = k. This 
form is not sufficiently general, because we have to keep all the four squark components 
independent variables; we can show numerically that K1K4 ^ K2K3. Moreover, the profile 
function h{r) introduced in the same paper turns out to be zero: as a consequence, the 
correct ansatz takes a much simpler form. The conclusions of Ref. |20j about the vortex 
moduli space although are not changed by these technical points. 




The ansatz for the gauge fields is: 



r 



r 



12 



gauge profile functions vanish and all the squark ones go to the value y/£/2. 
The boundary conditions at r — > are: 

f(r) = 2 + 0(r 2 ), f 3 (r) = (1 + cos a) + 0(r 2 ), g{r) = 1 + 0(r 3 ), (26) 

Kl (r)->(9(r 2 ), « 2 (r) -> O(r), « 3 (r) -> O(r), « 4 (r) -> O(l) . 
For the (2, 0) vortex we have: 

a = 0, 01 = Ki, 02 = «4, 

while for the (1,1) vortex (after a simple diagonalization): 

« = 7T, 0i = K 2 = 02 = «3- 

For the BPS vortex it is simpler to consider first order equations; but 
we are interested in understanding what is happening for 770,3 7^ 0. Thus 
we will write the equations in a form that can be easily generalized to a 
non-BPS setting. This will also give the possibility to check our equations 
and numerical results, just comparing the result for the tensions against the 
exact Bogomol'nyi bound; for completeness, we provide the first order BPS 
equations in Appendix B. The energy density due to the kinetic part of the 
gauge field is: 

, fj , sinW 
2r 2 eo 2 2r 2 e 3 2 2r 2 e3 
The part due to the kinetic energy of squark is 



S 9 - o„2 a 2 + Or2 a 2 + 0^2 o 2 ■ ( 27 ) 



rt 2 a I 12 1 l2 \ 1 • 2 a / 12 , I2\ , 

5 Q = COS 2 K + K 4 ) + S111 g ^ 2 + ^ ) + 

2 a /(l - cosa + /o + /s) 2 ^? . (1 - cosa - / + / 3 ) 2 ^1 



+ 2 1, S3 ^ + ^ 53 ^ 1 + 

. 2 « / (l+COSft-/o-/ 3 ) 2 ^ (l+COSQ + / -/ 3 ) 2 /t| \ 

+ Sm 2 V 4r 2 + 4r 2 1 



(1 — g) sin a 



4r 4 

(1 — g) sin 2 a 
2^ 



(cos 2 I (/t 2 + k\) + sin 2 I (k 2 + /t 2 )) - 

((1 + /o)kiK 3 + (1-/o)«2«4)- (28) 
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The part due to the potential reads: 




It is straightforward to write the Euler-Lagrange equations for this energy 
density, which are a system of seven second order equations, one for each 
profile function; for brevity we will not show them explicitly in the paper. We 
have solved this system numerically with the same method used in Sect. 3.2; 
in Fig. [3] is shown an example of the solution. The tension is found to be 
equal to T BPS = 4tt£ with an excellent precision for every a; this is a good 
numerical check for the solution obtained. 




2.5 5 7.5 10 12.5 15 17.5 2.5 5 7.5 10 12.5 15 17.5 



Figure 3: Vortex profile functions for a = tt/3. In the left panel there are k% (solid), 
K2 (long dashes), k 3 (short dashes), K4 (dots); in the right panel there are /o (solid), / 3 
(long dashes), g (short dashes). The following numerical values have been used: £ = 2, 
e - 1/4, e 3 = 1/2). 

An analytical check of the ansatz can also be found substituting Eqs. (j24|25p 
into the Euler-Lagrange equations (fl~2"]) . With this approach we find a sys- 
tem of the same seven second order equations with the following first order 
expression: 

K = elr 2 (K 3 K[ - kik' 3 + K 2 n i - k a k 2 ) - (1 - g)f z - (f 3 - cosa)g' = 0. (31) 
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This seems to be a paradox, because this is a system of eight differential 
equation with seven unknown functions. Actually, everything is consistent, 
because using the seven second order equations we can show the following 
property: 

f - ^ (32) 
dr r 

which shows that K is linear in r. From the boundary conditions of the 
profile functions, we find that the coefficient of this linear function has to be 
zero. This shows that our ansatz is consistent with the equations of motion. 

4.2 The non-BPS case: 7/0, 773 7^ 

For 770, r)3 7^ the (1, 1) and the (2, 0) vortices are still solutions to the equa- 
tions of motion; so these field configurations are extremal points of the energy 
which can be local minima or saddle points. For generic values of the pa- 
rameters, we find T^i) 7^ ^(2,0) j so the continuous moduli space interpolating 
between these two particular solutions disappears. For small values of r/o, 773 
we expect that the low energy physics of these solitons is described by an 
effective potential of the moduli space. In this section, we will estimate this 
potential numerically for generic values of a. 

A constraint on this potential comes from the BPS limit at r] ,r] 3 = 0. 
In this continuous family of degenerate solutions exists, with tension 

T = 47r£. If we insert these solutions into the energy density for 7/0, 773 7^ 0, the 
energy of these field configurations does not change. However, the solutions 
to the second order equations have energies which are less than this value. 
This sets an upper bound: 

T(a) < T BPS = 4tt£. (33) 
There is an obvious invariance of the equations: 

a — > —a. 

Indeed, if we expand around a = + <5ora: = 7r-|-5we find that the linear 
order in 5 is zero and that the first non-trivial correction to the tension is 
0(5 2 ). This shows that solutions with a = 0, n correspond to local minima 
or maxima of the tension. In order to find which of the two alternatives 
holds, an explicit calculation is needed. 
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In order to compute the potential of the vortex moduli space, we gener- 
alize the ansatz that we have used for the solutions in the BPS case, using 
the same expressions for the gauge fields, Q and the following expression for 
the adjoint fields: 

a-o = A (r), a 3 = A 3 (r), 
ai = (since)— Ai 2 (r), a 2 = (sina)— — Ai2(r), (34) 

where we have introduced the profile functions A , A 3 , Ai 2 , with the following 
boundary conditions: 

A (oo) = 0, A 3 (oo) = 0, A 12 (oo) = 0, (35) 

and the following r — > behaviour: 

AqocC(I), A 3 ocO(l), Ai 2 ocO(r). (36) 

This ansatz is suggested by the expression we get for these adjoint fields 
in the limit of large t]o,V3y where we can integrate these fields out (see Ap- 
pendix A). In the following we replace these expressions in the action and 
find second order equations for the profile functions for generic a. These field 
configurations at a ^ 0, 7r are not solutions to the full equations of motion, 
Eqs. ( TT2|) : they are just functional generalizations of the BPS solutions. We 
use these profiles as reasonable test functions to compute the effective moduli 
space potential. 




Figure 4: Vortex profiles for a = ir/3 and £ = 2, e = 1/2, e 3 = 1/4, r/ = rfy = 1. In 
the left panel there are K\ (solid), «2 (long dashes), K3 (short dashes), K4 (dots); in the 
middle panel /o (solid), f 3 (long dashes), g (short dashes); in the right panel Ao (solid), 
A3 (long dashes), A12 (short dashes). 

The kinetic energy of the adjoint scalars is: 

_ A /2 A 3 /2 + sin(a) 2 Ai 2 ' 2 sin(a) 2 ; 2 2 

^-71 + 12 + r 2 p 2 \{ L -9) A 3 + 

e e 3 ' e 3 
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+(/ 3 - cos(a)) 2 A 12 2 + 2 A 3 A 12 (cos a - / 3 ) {g - 1) } 



(37) 
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Figure 5: Vortex tension as function of a in the topological winding 2 sector; the tension 
of the BPS 2-vortex is normalized to Tbps = 2. In the left panel: £ = 2, eo = 1/2, 
e 3 = 1/2, r?o = 0.1, r?3 = 1. In the middle panel: £ = 2, e = 1/2, e 3 = 1/2, i] n = 1, 
773 = 0.1. In the right panel: £ = 2, e = 1/2, e 3 = 1/4, ij = V3 = 1< 

The potential term is: 



a 2 a 2 

cos(-) (/?! 2 + k 4 2 ) + sin(-) (K 2 2 + /«3 2 ) + 2r/oAo-^ ) + 

,a 2 



+ y ( cos (f )' ( K i 2 - K 4 2 ) + sin(|) (k 2 2 - k 3 2 ) + 2 773 A 3 ) + 

p 2 

+—sin(a) («i«3 - k 2 k 4 + 2r/ 3 Ai 2 ) +2sin(o;) A A i2 («i k 3 - k 2 k 4 ) + 



+ sinf 



J k 3 2 + cos(|) 2 k/^J ((Ao - A 3 ) 2 + sin(a) 2 A 12 2 ) + 
+ ( cos(|) 2 k, 2 + sin(|) 2 k 2 2 J ((Ao + A3) 2 + sin(a) 2 A 12 2 ) . (38) 



The energy is the sum of all pieces: 

£ = 2ttJ rdr(S g + 2S Q + S a + V). 



(39) 



The system of ten second order differential equations which we obtain is 
quite complicated, but can still be solved numerically. The qualitative plot 
of the profile functions is similar to the BPS case (see Fig. [3] and Fig. H|). 
The tension is a non-trivial function on the coordinate a, which gives us an 
effective potential of the moduli space. We solved this equations numerically 
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for different values of the couplings e ,e 3 , rj Q ,rj 3 and we have found three 
different regimes (see Fig. \5§. The tension can have a maximum at a = 
and a minimum at a = it; we also find the opposite situation in which there 
is a minimum at a = and a maximum at a = ir. The third alternative 
is that both a = 0, 7r are local minima of the tension, with one of them a 
metastable minimum. We never obtain a minimum at q ^ 0, ir. 



5 Vortex Interactions at Large Distance 
5.1 Vortex profiles at large distance 

At large distance r from the center of the vortex, the equations can be lin- 
earized and solved analytically. Let us introduce the following notation: 



<h = VU^ + S<h., <h = Vt/2 + 8<h, v = (50i,<50 2 ,A o ,A 3 ). 
The following linear differential equations can be written: 



v"{r) + ^ - Wv{r) 



0. 



(40) 



(41) 



where the matrix W is given by: 

2 2 



w 



2 



2 



f e §^3 



2^e| + ^ ) 



(42) 



The eigenvalues of the matrix W are in direct correspondence with some of 
the scalar spectrum of the theory (see Eq. (ITU]) ): 



(43) 



^3,4 = M| 1)T2 = £e 2 + e^ 3 2 ± ^/2^3 2 e 6 + e^f . 
The corresponding eigenvectors are: 



Ul,2 



e W ± v / 2^o 6 r ?0 2 + ep^o 4 -e W ± v^ep^o 2 + ep^o 4 



2^eo 4 r ?0 



2^2£e 4 7 ? o 



.0 
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- _ ( -e 3 W ± y / 2ee 3 6 r?3 2 + e 3 N? e 3 W =j= VggV + ^8^4 \ 
V3 ' 4 ^ 2^/2^3 ' 2^63^3 ' U '7' 

Note that v k are also eigenvectors of the mass matrix defined in Eq. (Q. 

The solutions to these equations which are zero at infinity are given by 
the modified Bessel function: 



b k v k K (yftF k r), (44) 



v{r) = 

k=l,-,4 

where h k are appropriate constants which can be found solving the complete 
differential equation also at small r. For large x we can use: 

K {x) « y^e-. (45) 
The asymptotic solutions for the scalar profiles read: 



fc=l,-,4 

The large r equations for / 3 and /o are: 

/s f -^-e^/o = o, ^-^-^ /3 o. (47) 

This leads to the following asymptotic expression in terms of Bessel functions: 

/o, 3 = c^rK^e^r) cc y^e~^^ r , (48) 

where cq, c 3 are constants which should be determined by the original 2nd 
order differential equations. Note that there is the identity: K\{r) = —K' Q (r). 



5.2 Static vortex potential 

The next step is to reproduce the vortex asymptotic interactions in the ef- 
fective linear theory by coupling the low-energy degrees of freedom to an 
effective scalar density p and an effective vector current j M . We shall extend 
an approach used in Refs. [3TJ [32] . 

First of all, we have to discuss the bosonic particle spectrum of the theory. 
There is a massive U(l) vector and a massive SU(2) vector; then in principle 



19 



there are 12 complex scalar fields (Q, Q, A, A k ). For the vortex solution we 
have used the ansatz Q = Q\ so in order to discuss the vortex interactions 
we can neglect the modes that break this condition^]. There are 16 real fields, 
4 of which are eaten by the Higgs mechanism; finally we have 12 physical 
scalars. We have already calculated the masses of these particles in Sect. 2.2; 
at low energy and at low coupling we can write a free theory which describes 
the infrared physics: 

Zfree = W + |^ + ^| W + ^K A ^ + ^ 

+^) 2 + \{d^f + ^Sf + ^f{T k f. 

This effective Lagrangian contains three real scalar fields, $=0,1,2, which are 
SU(2) singlets and three real scalar fields, T k =0 1 2 , which are SU(2) triplets. 
These scalars correspond to the appropriate eigenvectors of the mass ma- 
trix in Eq. (jUj). The index k is an SU(2) triplet index; this SU(2) group 
corresponds to the SU{2) c +f in the full theory. 

In order to include external sources (=vortices) in this effective Lagrangian, 
we need the following effective terms: 

Csource = PsA + P^T* + + j^A*". (50) 

The corresponding wave equations are: 

(D + M^S^psi, (□ + M 2 Tl )T k = p k Tl , (51) 

(□ + e 2 0^° M = J M , (□ + e 3 2 0^ = j k "- 

On the other hand, for the (1, 0) vortex with orientation n k , we have the 
following asymptotic profiles, converted into the singular real Q gauge: 

S = 0, S 1 = hK (M Sl r) S 2 = b 2 K (M S2 r), (52) 

T fc = 0, T k = b 3 n k K (M T1 r) T k = b 4 n k K (M T2 r) , 

A = -co(z A VK (e ^r)), A k = -c 3 n k (z A VK (e 3 ^r)), 

9 If we wish to include these extra modes in the low energy theory, we need only to 
promote the real fields Si,S2, T^, T k in Eq. (|49|) to complex fields. 
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where V is only the ordinary gradient and not the covariant derivative as in 
the other sections. We need the following mathematical identity for the 2 + 1 
dimensional Laplacian of Kq in term of Dirac's 8 function: 

(-A + M 2 )K (Mr) = 2n5(r). (53) 

The following expressions are found for the scalar densities corresponding to 
a vortex placed at the position x and having orientation n k : 

pso = 0, p S i = 2nb 1 6(x), p S 2 = 2nb 2 S(x), (54) 

Pro = 0, pri = 2irb 3 n k 5(x), p T2 = 2nb 4 n k 8(x). 
In a similar way we obtain the following expressions for the currents: 

J = -2irc z A V5(£), J k = -2nc 3 n k z A V£(x). (55) 

Using these expressions, it is straightforward to compute the static inter- 
vortex potential between two vortices with orientations n\ and n 2 at distance 
R: 

U = 2n (c 2 K (e ^R) - b\K Q (M sl R) - b 2 2 K (M S2 R)+ (56) 



+ (n 1 ■ n 2 )(c 2 3 K (e 3 ^R) - b 2 3 K (M T1 R) - b 2 4 K (M T2 R))) . 



In the BPS case (r/o = i]?> = 0) this potential is exactly zero, because we have 
M S i = M S2 = eoV?, M T1 = M T2 = e 3 ^ and c 2 = bj + b 2 , c 2 = b 2 + b\. 



If i] 3 , r/o 7^ 0, we find that at large distance the particle with lowest mass 
is the one which dominates the interaction. We have always the following 
inequalities: 

M S2 < eov 7 ^ < M S i, M T2 < e 3 ^l < M Tl . (57) 
Thus if M S2 < M T2 then we have: 

U » -2nb 2 2 K (M S2 R) » -2nb 2 2 ^ ^-^e~ Ms2R , (58) 

which gives an always attractive force. On the other hand, if Mt 2 < Ms 2 : 

U » -27r^X (M T2J R)(n 1 • n 2 ) » -2nb 2 I— e^Vi • n 2 ), (59) 

y ZlVlj' 2 rC 
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which gives attraction for n% = n 2 and repulsion for n% = —n 2 . 

A very peculiar thing happens for e$ = and J)o = % ^ 0. For these 
fine-tuned values of the couplings Ms2 = Mt2 = M 2 and b 2 = 64, so the 
effective vortex potential has the form: 



U « -2nb 2 2 K (M 2 R)(l + n t ■ n 2 ) « -2nb 2 2y J ^-^ e - M * R (l + n x ■ n 2 ), (60) 

which gives a flat potential for Hi = —n 2 . This is consistent with the fact 
that in this limit the (1, 0) and the (0, 1) vortices do not interact because they 
are completely decoupled (see the argument below Eq. f[2"2l ). This behaviour 
is similar to the one found in Ref. [25] for global non-Abelian vortices. 



6 Effective worldsheet theory 
6.1 Single vortex 

It is useful in the following to use the singular gauge in which the squarks 
fields at r — > 00 tend to a fixed VEV and do not wind. In this gauge, the 
ansatz ([IE]) for the single vortex reduces to 

A° = »/ 0l A{ = '-^JW. a = A , a k = n%, (61) 

Q = ^±^1 + " Vn fc . 
v 2 2 

We will assume that the orientational coordinates n are functions of the 

string worldsheet coordinates, n becomes a field of a 1 + 1 dimensional sigma 

model. This effective theory has no potential due to the fact that the n k 

parameterize some zero modes; in the following we will compute the kinetic 

term. For the gauge field components A ,3, we will use the same ansatz as 

used in Refs. [6j [T6] : 

A k = -^(r a e abc n b d k n c )p(r), k = 0,3. (62) 

The field strength components F k i with k = 0, 3 and i = 1, 2 are not zero 
any more: 

F kl = ~ d k n a r a e l3 ^ / 3 [1 - p(r)] + ^( T a e abc n b d k n c )^^-p(r). (63) 
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Substituting this expression into the kinetic term for the gauge field and 
for the squark fields, we obtain a simple generalization of the BPS case 
discussed in Refs. [Bl [TBI E3] : 

£1+1 = £ j dtdz ( d . n ky } ( 64) 

where: 

P=%frdr \f + 7 Pf + At(l - p) 2 + (65) 

+e^{(0? + ^)^ + (l-p)(0i-0 2 ) 2 }}- 

We have to solve the Euler-Lagrange equations for p(r), with p(0) = 1 and 
p(r — > oo) = 0. In the BPS case, where A3 is trivially 0, we can show 
from the equations of motion [16j that p = 1 — <\>\j<\>2 and that (3 = 2%/e\ 
(see also [331 131]); in the general case r] 3 ,r] 7^ there is not such powerful 
analytical result, here we have to solve the equations for p numerically and 
then calculate (3. 

In the BPS case we have additional fermionic zero modes, associated with 
the unbroken supercharges; at small r/0,^3 7^ these modes should be still 
present, but they will not be described anymore by the fermionic sector of a 
super symmetric effective theory in 1 + 1 dimensions. We will not discuss this 
aspect in this paper and we will leave it as a problem for further investigation. 

The color-flavor modes of the (2,0) vortex are very similar to the (1,0) 
ones: both the vortices have a CP 1 moduli space and the value of (3 can be 
determined using Eq. (|65|) . For the (1, 1) vortex,nevertheless, these modes are 
just trivial because all the profile functions are proportional to the identity 
matrix. 

6.2 Two well separated vortices 

A proper description of the system has to take into account also the quantum 
aspects of the sigma model physics. Let us consider two vortices with internal 
orientations ni,^. The relative distance between them can be promoted to 
a complex field R; the global position of the system, on the other hand, 
decouples from the other the degrees of freedom. If the distance of the 
two vortices is large (\R\ — > 00), we expect that the effective worldvolume 
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description of the bosonic degrees of freedom is: 



(66) 

where: 



where T is the tension of a single vortex. This description is good only for 
large values of the VEV of the field R; at R — the internal degrees of 
freedom are no longer described by CP 1 x CP 1 , but by a space with topology 
(see Refs. [20| l2T]). Moreover, the expression used for the potential 
is good only for large vortex separation. 

If we keep the VEV of R fixed (which physically corresponds to keep the 
distance of the two vortices fixed with some external device), the effective 
description is given by two CP 1 sigma models with a small interaction term 
of the form cn\ ■ n 2 . 



7 Conclusion and Discussion 

For Abelian type I superconductors, the force between two vortices with the 
same winding number is always attractive. This is true at large and at small 
distances, as shown by numerical calculations in Ref. [21]. In the model 
discussed in this paper, for r/ , r/ 3 > 0, the masses of some of the scalars fields 
are always found to be less than the mass of the corresponding vector boson. 
In this sense we can think the system as a generalization of the Abelian type I 
superconductor. However, here is an important difference: the force between 
two vortices is not always attractive; there is a non-trivial dependence on the 
coupling, the relative internal orientation and the distance. 

In this paper, we have studied the problem in two different limits: large 
vortices separation and coincident vortices. For large separations we have 
computed the leading potential analytically; the behavior at large distance 
is dominated by the particle with the lowest mass M low . There are two main 
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alternatives, which hold for different values of the couplings: 



U(R) oc 



2M S2 R 



1 P -M S2 R 



for Mi 



low 



M S2} Type I 



(68) 



-{nx ■ n 2 




low 



M T2 , Type I* 



where Ms 2 , M?2 are the masses of the scalars in Eq. fflOl) . In order to 
distinguish these regimes, we call them Type I and Type I*; for Type I* 
vortices the sign of the asymptotic force depend on Hi ■ H 2 . For the fine- 
tuned values e = e 3 and i] = r/ 3 7^ 0, the relation M S2 = M T2 = M 2 holds, 
and the effective vortex potential has the form: 



which gives a flat potential for Hi = —n 2 . 

For coincident vortices we have found two stationary solutions of the 
equations of motion, the (1,1) and the (2,0) vortices, and we computed 
their tensions numerically. The results are shown in Fig. [2j both the cases 
Ti t i > T 2) o or T 2i o > T^i are possible for different values of the coupling. The 
moduli space interpolating between these solutions at r/ = r/ 3 = disappears 
for non-zero values of one of these parameters (see Fig. [5]). 

It is interesting to match the data of the two complementary approaches. 
Let us for simplicity consider the case of parallel (Hi = n 2 ) and anti-parallel 
(Hi = — n 2 ) vortices. In the case of parallel vortices at large separation 
distance, the force is always attractive; also from numerical calculations we 
find T 2 < 2T 10 for all the values of the coupling that we have analyzed. We 
have not made the calculation for arbitrary distances, but we think that the 
above are a good evidence for the fact that the force between two parallel 
vortices is always attractive in our model. 

For anti-parallel vortices, on the other hand, the situation is more com- 
plicated. At large distances there is attraction if Mti > Ms 2 and repulsion 
if M T2 < M S2 . For the choice e = e 3 and i] = r/ 3 , the relation M T2 = M S2 
holds; as we already noticed in Sect. 3, for these particular values the two 
diagonal U(l) factors decouple: there is no net classical force between the 
(1, 0) and the (0, 1) force for arbitrary distance (this configuration is not sta- 
ble, because if we allow Hi, n 2 to vary, we have that (2, 0) vortex has a lower 
energy). Indeed, if we keep e = e 3 , we obtain M T2 > M S2 , T11 < 2T 10 




(69) 
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Figure 6: Qualitative plot of the vortex potential as function of the vortex distance for 
fti = —ri2 and eo = e$. For 770 > 773 we have attraction (left); for 770 < 773 we have repulsion 
(right). With the choice rjo = 773 there is no net classical force. 

for 770 > 773 and Mti < Mg 2 , ^1,1 > 2Ti )0 for rjo < r/ 3 . This is a good evi- 
dence that for rjo > 773 we have an attractive force and for 770 < 773 we have a 
repulsive one (see Fig. (jSJ)). 




Figure 7: For eo 7^ e% the qualitative plot of the vortex potential as function of the 
vortex distance for n\ = —712 can have maxima or minima for non-zero vortex separation. 

If we relax the condition eo = e 3 , there are situations in which at large 
distance there is attraction (because M T2 > M S2 ) and also we get > 
2T 10 , as is shown in Fig. [7| to the left: this means that there is a critical 
distance, in which there is a minimum of the inter-vortex potential for anti- 
parallel vortices (if we allow the vortex orientation to flip, probably it will 
not be a minimum any more, because we still have T 2 fl < 2T\$). An example 
of this situation can be obtained with the couplings £ = 2, eo = 1/4, = 
1/2, rjo = 0,773 = 4. Moreover we can obtain Mt2 < M$ 2 and T\ t \ < 2Ti i0 , 
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which means that there is a critical distance at which there is a maximum 
of the inter- vortex potential (see Fig. [7] on the right). An example of this 
situation can be obtained with the couplings £ = 2, e = 1/2, e% = 1/4, rjo = 
3.5,773 = 2. 

It is interesting that at R — there are two different regimes, which 
depend on the values of the coupling, with very different properties. The 
physics of the (2, 0) vortex is described by a bosonic CP 1 sigma model; the 
(1,1) vortex on the other hand is an Abelian vortex with no internal de- 
grees of freedom. For some values of the coupling we have an evidence that 
both vortices are local minima of the tension (see Fig. [5]): one of the two is 
metastable (indeed for some fine tuned values of 770,3 we have that both the 
vortices have the same tension). 

A model with metastable vortices at weak coupling has already been 
studied in Ref. [35J. This behavior is reminiscent of SU(N) Yang-Mills, 
where for each topological n-ality we can have different string tensions for 
each representation of the Wilson Loop. In each topological sector there 
is just one stable string, corresponding to the antisymmetric representation; 
there are evidences that the strings with other representations are metastable 
strings, at least in the large N limit (see Ref. [36] for a discussion). 

In this paper we have been interested in non- Abelian non-BPS vortices in 
an M = 1 supersymmetric model. In such a restricted model we have found 
a physics similar to Abelian type I superconductors. In a companion paper, 
we will discuss a simpler non-supersymmetric model in which we can have 
both type I and type II non- Abelian superconductivity. 
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A Large 770,773 limit 

If we take the limit e ^o \/£ we can integrate out the superfield a from 
the superpotential: 

a = — . 70) 

Similarity way in the limit e 3 i] 3 » y/£ we can integrate out the superfield 



a k : 



j. = -Tw^y (71) 

2>(3 



The effective superpotential is: 

Tr(Qr k Q)TT(Qr k Q) + (£-Tr(gg)) 2 j 



4^3 4?7 



The potential is: 



F = l^gV^-Tr^r^Vl^gt^-Tr^))^ 



|e-TrgQ| 2 , |TrQr fc Q| 2 



+ 1 „ o 1 +- — ^ + 



8^ 8^ : 



Tr(gtg) + Tr(QQt)) 



+ 



i6 klm Tr(QW k &)Tr(Qr l Q) 



(Tr(Qt r -Q) + Tr(gr m gt) j _. 



(g - Trgtgt) Tr (gr fc g) + g - TrgQ)Tr(QVfcgt) 

(Tr(gtr c g) + Tr(gr c gt)) . (73) 

Note that in this low energy action there is another vacuum at Q = Q = 0. 

The equations of the (p, k) vortices are a bit simpler, because we can 
integrate out the adjoint fields a,a k and so we need less profile functions. 
The energy is: 
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, (# + 02)(/ o 2 + /f) + 2/ 3 /o(^ - g) , (0j + - # - 0IO 2 , 



2r 2 Ar]l 
(0 2 -0l) 2 (0 2 + 0l) (e-0?-0l)(0 2 -^) 2 



4ry| 2?7o?73 

The corresponding Euler-Lagrange equations are very similar to Eqs. (|20|) . 

Numerical calculations can be performed for the profile functions and the 
tension. In Fig. [H] there is a comparison between the tension calculated in 
the full theory and in the large 770,3 approach. For small rjj the correction in 
the tension from the BPS case is quadratic in rjj, as discussed in the Abelian 
case in Ref. 



Figure 8: Tensions for the (1,0) for 773 = rjo = r). The numerical values eo = 1/2, 
e3 = 1/4, £ = 2, < 77 < 10 are used. The black dots give the result of the calculation in 
the full theory; the blue one give the result in the large rjj effective theory. For rjj = the 
vortex is of course BPS. 



The asymptotic profiles are also simpler in this limit. Let us define: 

S = (501 + (502, d = 501 — 502, 

we the find two linear differential equation: 

8 » + 1 _ 5 = 0, d" + - - %d = 0. 
r r/Q r V3 

So the asymptotic solutions for the field profiles are: 

-(V2£A) 3 )r 



(75) 



(76) 



'1,2 



«0 



(77) 



where s , d are constants analogous to b^. 
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B BPS equations 

In terms of fields, the BPS equations for two coincident vortices for rj Q , i] 3 = 
read: 

F12 + e 2 3 Tr(Q VQ) = , F° 2 + e\ (Tr(Q*Q) -0 = 0, (78) 
(Vi + iV 2 )Q = 0. 

In term of profiles functions the following system of seven first order equations 
holds: 



fo 
r 



e ° { ( c ° s f ) 2 ^ + ^ + ( sin f ) 2 ^ + ^ ~ 4 ' (79) 

f = e^{(cos|) 2 K-^)+(smf) 2 (^ ' 2 



5 2 



r 



/ 0-1 . 2 /«\ , 1 - cos(a) + f + fo 
Kl = — sm (- j % + ^ «i, 

, - 1 2 f a \ 1 + cos(a) - f - h 

K 2 = —COS ^-J« 4 ^ K 2 , 



/ 5-1 2 ( a \ , 1 + cos(a) + f - fo 

K 3 = — — COS (2)^ + ^ K 3, 

/ 5-1 . 2 ( a \ 1 - cos(a) - f + fo 

k a = -— sm [- ) K 2 K A . 

In our numerical analysis of Sect. 6.1 we used the second order Euler-Lagrange 
equations; at the end we used these first order equations as a check of our cal- 
culation. We have found an excellent agreement between the two approaches. 



References 

[1] G. 't Hooft, Nucl. Phys. B 190 (1981) 455. S. Mandelstam, Phys. Lett. 
53B (1975) 476; Phys. Rep. 23C (1976) 245. 

[2] A. Abrikosov, Sov. Phys. JETP 32 1442 (1957), H. Nielsen and P. Ole- 
sen, Nucl. Phys. B 61 45 (1973). 



30 



[3] N. Seiberg and E. Witten, Nucl. Phys. B 426 (1994) 19, 



hep-th/9407087 



[4] M. R. Douglas and S. H. Shenker, "Dynamics of SU(N) supersymmetric 
gauge theory," Nucl. Phys. B 447 (1995) 271 | arXiv:hep-th/9503163] ; 
A. Hanany, M. J. Strassler and A. Zaffaroni, "Confinement and strings 
in MQCD," Nucl. Phys. B 513 (1998) 87 |arXiv:hep-th/9707244| . 



[5] A. Hanany and D. Tong, "Vortices, instantons and branes," JHEP 0307 
(2003) 037 |arXiv:hep-th/0306150| . 

[6] R. Auzzi, S. Bolognesi, J. Evslin, K. Konishi and A. Yung, "Nonabelian 
superconductors: Vortices and confinement in N = 2 SQCD," Nucl. 
Phys. B 673 (2003) 187 |arXiv:hep-th/0307287| . 



[7] D. Tong, "TASI lectures on solitons," arXiv:hep-th/0 509216; M. Eto, 
Y. Isozumi, M. Nitta, K. Ohashi and N. Sakai, "Solitons in the Higgs 
phase: The moduli matrix approach," J. Phys. A 39 (2006) R315 
|arXiv:hep-th /0602170]. K. Konishi, "The magnetic monopoles seventy- 
five years later," arXiv:hep-th/0702102; M. Shifman and A. Yung, "Su- 
persymmetric Solitons and How They Help Us Understand Non-Abelian 
Gauge |arXiv:hep-th/0703267| 



[8] M. Eto, K. Hashimoto and S. Terashima, "Solitons in super- 
symmety breaking meta-stable vacua," JHEP 0703, 061 (2007) 
|arXiv:hep-th/0610042|. M. Shifman and A. Yung, "Confinement in 
N=l SQCD: One Step Beyond Seiberg's Duality," larXiv: 0705 .38111 
[hep-th]; A. Gorsky, M. Shifman and A. Yung, "N = 1 supersymmet- 
ric quantum chromodynamics: How confined non-Abelian monopoles 
emerge from quark condensation," Phys. Rev. D 75 (2007) 065032 
[arXiv:hep-th/0701040| ; M. Eto, K. Hashimoto and S. Terashima, 
larXiv:0706.2005l [hep-th] . 

[9] M. Eto et ai, "Non-Abelian duality from vortex moduli: a dual model 
of color-confinement," arXiv : hep-th /0611313. 

[10] L. Ferretti, S. B. Gudnason and K. Konishi, "Non-Abelian vortices and 
monopoles in SO(N) theories," larXiv:0706.3 854 [hep-th]. 

[11] A. Gorsky and V. Zakharov, "Magnetic strings in Lattice QCD 
as Nonabelian Vortices," larXiv:0707.1284l [hep-th]; A. Gorsky and 



31 



V. Mikhailov, "Nonabelian strings in a dense matter," arXi v:0707.2 304 
[hep-th] . 

M. Shifman and A. Yung, "Non-abelian flux tubes in SQCD: Super- 
sizing world-sheet supersymmetry," Phys. Rev. D 72 (2005) 085017 
[arXiv:hep-th/0501211 . 



M. Edalati and D. Tong, "Heterotic vortex strings," JHEP 0705 (2007) 



005 |arXiv:hep-th7 0703045|; D. Tong, |arXiv:hep-th/0 703235, 



A. Gorsky, M. Shifman and A. Yung, "Non-Abelian Meissner effect in 
Yang-Mills theories at weak coupling," Phys. Rev. D 71 (2005) 045010 
|arXiv:hep-th/0412082| . 

D. Ton g, "Monopoles in the Higgs phase," Phys. Rev. D 69 (2004) 
065003 |arXiv:hep-th/0307302] . 

M. Shifman and A. Yung, "Non-Abelian string junctions as confined 
monopoles," Phys. Rev. D 70 (2004) 045004 |arXiv:hep-th/0403149] . 

A. Hanany and D. Tong, "Vortex strings and four- dimensional gauge 
dynamics," JHEP 0404 (2004) 066 |arXiv:hep-th/0403158| . 

M. Eto, Y. Isozumi, M. Nitta, K. Ohashi and N. Sakai, "Moduli 
space of non-Abelian vortices," Phys. Rev. Lett. 96 (2006) 161601 
arXiv:hep-th/0511088|. 



K. Hashimoto and D. Tong, "Reconnection of non-abelian cosmic 
strings," JCAP 0509 (2005) 004 |arXiv:hep-th/0506022| . 

R. Auzzi, M. Shifman and A. Yung, "Composite non-Abelian 
flux tubes in N = 2 SQCD," Phys. Rev. D 73 (2006) 105012 
|arXiv:hep-th/0511150l . 

M. Eto, K. Konishi, G. Marmorini, M. Nitta, K. Ohashi, W. Vinci and 
N. Yokoi, "Non-Abelian vortices of higher winding numbers," Phys. Rev. 
D 74 (2006) 065021 |arXiv:hep-th/0607070| . 

M. Eto, K. Hashimoto, G. Marmorini, M. Nitta, K. Ohashi and 
W. Vinci, "Universal reconnection of non-Abelian cosmic strings," Phys. 
Rev. Lett. 98, 091602 (2007) larXiv:hep-th/0609214| . 



32 



[23] E. Bogomolny, The stability of classical solutions, Sov. J. Nucl. Phys. 
24, 449 (1976). [Reprinted in Solitons and Particles, Eds. C. Rebbi and 
G. Soliani (World Scientific, Singapore, 1984), p. 389]. 

[24] L. Jacobs and C. Rebbi, "Interaction Energy Of Superconducting Vor- 
tices," Phys. Rev. B 19 (1979) 4486. 

[25] E. Nakano, M. Nitta and T. Matsuura, "Interactions of Non-Abelian 
Global Strings," larXiv:0708.4092l [hep-ph]; "Non-Abelian Strings in 
High Density QCD: Zero Modes and Interactions," larXiv:0708.4096l 
[hep-ph] . 

[26] S. C. Davis, A. C. Dav is and M. Trodden, Phys. Lett. B 405 (1997) 257 
[arXiv:hep-ph/9702360|; A. Gorsky and M. A. Shifman, Phys. Rev. D 
61 (2000) 085001 |a rXiv:hep-th/9909015] . 

[27] A. Achucarro and T. Va chaspati, Phys. Rept. 327, 347 (2000) [Phys. 
Rept. 327, 427 (2000)] |arXiv:hep-ph/9904229| . M. Eto et al, "On 
the moduli space of semilocal strings and lumps," larXiv:07 04.2218 
[hep-th]. M. Shifman and A. Yung, "Non-Abelian semilocal strings 
in N = 2 supersymmetric QCD," Phys. Rev. D 73, 125012 (2006) 
|arXiv:hep-th/0603134] . 

[28] A. I. Vainshtein and A. Yung, "Type I superconductivity upon monopole 
condensation in Seiberg-Witten theory," Nucl. Phys. B 614 (2001) 3 
|arXiv:hep-th/001225"0l . 

[29] R. Auzzi, S. Bolognesi and J. Evslin, "Monopoles can be confined by 0, 1 
or 2 vortices," JHEP 0502 (2005) 046 |arXiv:hep-th/0411074| ; S. Bolog- 
nesi, "The holomorphic tension of vortices," JHEP 0501 (2005) 044 
JarXiv:hep-th /041 1075| ; S. Bolognesi and J. Evslin, "Stable vs unstable 
vortices in SQCD," JHEP 0603 (2006) 023 |arXiv : hep-th/ 05 06 174| . 

[30] M. Shifman and A. Yung, "Localization of non-Abelian gauge fields on 
domain walls at weak coupling (D-brane prototypes II)," Phys. Rev. D 
70 (2004) 025013 |arXiv:hep-th/0312257] . 

[31] J. M. Speight, "Static intervortex forces," Phys. Rev. D 55 (1997) 3830 
|arXiv:hep-th/9603155l . 



33 



[32] A. Marshakov and A. Yung, "Non-Abelian confinement via Abelian flux 
tubes in softly broken N = 2 SUSY QCD," Nucl. Phys. B 647, 3 (2002) 
|arXiv:hep-th/0202172] . 

[33] M. Eto, Y. Isozumi, M. Nitta, K. Ohashi and N. Sakai, "Manifestly 
super symmetric effective Lagrangians on BPS solitons," Phys. Rev. D 
73, 125008 (2006) |arXiv:hep-th/0602289l . 

[34] M. Eto, Y. Isozumi, M. Nitta, K. Ohashi and N. Sakai, "Instantons in the 
Higgs phase," Phys. Rev. D 72, 025011 (2005) | arXiv:hep-th/0412048] . 

[35] M. Shifman and A. Yung, "Metastable strings in Abelian Higgs models 
embedded in non-Abelian theories: Calculating the decay rate," Phys. 
Rev. D 66 (2002) 045012 |arXiv:hep-th/0205025| . 

[36] A. Armoni and M. Shifman, "Remarks on stable and quasi-stable k- 
strings at large N," Nucl. Phys. B 671 (2003) 67 | arXiv:hep-th/0307 020|. 

[37] X. r. Hou, "Abrikosov string in N = 2 supersymmetric QED," Phys. 
Rev. D 63, 045015 (2001) [arXiv:hep-th/0005119] . 



34 



